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Section  1.  INTRODUCTION 


The  optimum  linear  filter  developed  by  Kalman  and  others  is  a  dynamic 
means  of  data  processing  which  is  useful  in  estimating  the  state  of  linear 
systems  from  noise  measurements  (Ref.  1).  Among  the  many  problems  which 
can  be  formulated  in  a  manner  compatible  with  the  Kalman  Filter  is  the 
transfer  alignment  of  one  inertial  guidance  platform  (slave)  to  another 
(master).  If  it  is  assumed  that  alignment  to  within  a  few  degrees  is 
possible  by  a  means  such  as  gimbal  angle  matching,  the  differential 
equations  which  describe  state  variable  differences  between  the  two  sys¬ 
tems  are  basically  linear.  The  differences  which  become  the  state  var¬ 
iables  of  the  alignment  problem,  are  random  variables  as  are  the  inertial 
sensor  errors  which  cause  them  to  change  with  time.  Certain  of  these 
alignment  state  variables  are  readily  available  to  use  as  inputs  to  the 
Kalman  Filter — that  is,  they  are  the  measurements.  The  misalignment 
angles  between  the  two  systems  can  then  be  estimated  and  removed  with  an 
accuracy  greater  than  that  obtainable  through  practical  mechanical  trans¬ 
fer  alignment  techniques. 

This  report  develops  base-line  data  for  use.  in  evaluating  proposed 
equipment  and  estimation  schemes  for  the  transfer  alignment  problem.  The 
principal  quantities  of  interest  are  the  Kalman  Filter  errors  in  estimates 
of  the  misalignment  angles.  It  is  assumed  that  correction  of  the  mis¬ 
alignment  can  be  easily  carried  out  with  a  greater  accuracy  than  the 
estimation,  therefore  estimation  accuracy  is  equivalent  to  alignment 
accuracy. 

Section  2  provides  a  derivation  of  the  alignment  system  equations, 
including  all  relations  necessary  to  analyze  the  Kalman  Filter  accuracy. 
Section  3  identifies  those  alignment  system  parameters  which  are  most 
important  and  Section  4  provides  trade-off  relations  between  them  and 
alignment  time  and  accuracy.  More  complex  schemes  aie  briefly  discussed 
in  Section  5  as  the  basis  for  future  investigation. 
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Section  2.  DERIVA‘i'*u.  Cl  SYSTEM  AND  MEASUREMENT  EQUATIONS 


BASIC  ERROR  DYNAMICS 

The  linear  equations  describing  the  dynamics  of  misalignment  angles 
and  other  Important  parameters  are  derived  in  this  chapter  for  the  transfer 
alignment  between  two  nearly  aligned  inertial  platforms.  The  systems  are 
assumed  to  be  aligned  to  within  a  few  degrees  before  the  start  of  the 
Kalman  Filter  operation.  Glmbal  angle  matching  techniques  are  available 
to  perform  this  coarse  alignment  and  the  repldual  attitude  errors  are 
assumed  to  result  from  flexure  of  the  structure  connecting  the  two  plat¬ 
forms  and  from  mounting  errors.  The  measurements  permitted  are  differences 
between  horizontal  velocity  indications  by  the  two  systems  operating 
independently  during  filtering.  It  is  assumed  that  velocity  corrections 
are  applied  to  account  for  the  fact  that  the  two  systems  may  be  physically 
separated  by  a  rotating  rigid  structure.  The  slave  platform  is  assumed 
to  be  provided  with  initial  values  of  velocity  and  position  from  the 
master,  and  the  earth  rate  platform  commands  applied  to  the  two  systems 
are  identical. 

It  can  be  seen  from  Eq.  42,  45,  and  50  of  Ref.  1  that  ideal  accel¬ 
erometers  measure  the  quantity 


(1) 


The  vector  of  master  system  accelerometer  outputs,  identified  by  the 
subscript  M,  is 

master  ace  terometer  outputs  ■(  —  )+■  VM  (2) 

\»/M  M 

where  is  the  vector  of  master  system  accelerometer  errors.  Because 
the  angular  misorientation  between  the  slave  accelerometer  input  axes 
and  the  corresponding  master  axes  is  small  it  can  be  represented  by 
a  vector  of  rotation,  The  components  of  £x*  £y.  and  £r,  represent 
th»  angular  rotation  about  the  master  axis  designated  by  the  e'^cripr 
which  is  necessary  to  bring  the  master  platform  Into  alignment  with  the 
slave  platfDrn.  Since  the  two  systems  are  closely  aligned,  the  negative 
of  these  componeuts  describes  the  rotation  about  the  designated  slave 
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axis  necessary  Co  bring  Che  slave  system  inco  alignment  with  Che  master 
(Che  order  of  roCation  is  unimportanc) .  NeglecCing  effects  due  to 
separation  between  the  two  systems,  the  vector  of  slave  accelerometer 
outputs  would  be 


slave  accelerometer  outputs 


(3) 


where  Vs  is  the  vector  of  slave  accelerometer  errors. 

It  is  assumed  that  the  differences  between  slave  and  master  system 
velocity  indications  are  corrected  for  the  relative  velocity  which  results 
from  their  physical  separation  on  a  rotating  rigid  body.  This  correction 
can  be  computed  from  knowledge  of  the  aircraft  angular  velocity  and  the 
relative  position  of  the  two  platforms  with  respect  to  the  aircraft  center 
of  mass.  Any  errors  in  this  calculation,  as  well  as  structural  flexure 
effects,  will  be  viewed  as  measurement  errors  in  the  velocity  difference. 
As  a  consequence  of  the  velocity  correction  applied,  the  difference 
between  Eq.  2  and  3  can  be  taken  as  the  derivative  of  the  noise-free 
velocity  difference  vector,  Av 


t 


Av  - 


(A) 


If  the  ideal  coordinate  axes  are  aligned  north  (x) ,  east  (y),  and  down  (z) 
and  if  the  master  system  is  closely  aligned  to  them,  we  get  (Eq.  1) 
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where  ax,  ay,  and  az  are  accelerations  relative  to  the  earth,  vx,  vy,  and 
v2  are  velocities  relative  to  the  earth  and  ftx  and  ftz  are  components  of 
the  earth's  rotation  vector.  The  magnitude  of  the  terms  arising  from 
Coriolis  acceleration  is  very  small.  For  example,  if  the  aircraft  has 


3 


NWC  TP  4653 


a  velocity  of  1,000  ft/sec,  the  magnitude  of  QxVy  cannot  exceed 
7.3  *  10“  2  ft/sec2.  Of  course  this  is  such  less  than  the  magnitude  of  g. 
Similarly,  ax  and  ay  will  be  much  larger  than  the  corresponding  x-axis  and 
y-axis  Coriolis  terms  unless  a  very  straight  flight  path  is  maintained. 
Consequently,  though  it  is  not  necessary  to  permit  the  use  of  the  Kalman 
Filter  in  transfer  alignment,  the  Coriolis  acceleration  terms  in  Eq.  5 
are  neglected  and  4  becomes 


a 


x 


Av  -  -£  * 


a 

.  z 


g 


(6) 


The  components  of  Av  are  given  by 

‘  (\  '  \)  *  *y  '  K  '  *)  S 
4ir  •  (7rs  -  \)  *  (*.  -  *>  -  *» 

•  (%  -  \)  \  °> 

The  two  platforms  experience  angular  rates  with  respect  to  inertial 
spjee,  expressed  in  their  own  coordinate  frames,  hy 


(“*)„  •  5  ’  '  S 

(“s)s  ■  5  *  5:  *  's 

(8) 

It  ie  assumed  that  the  earth  rate  rotational  commands  are  the  same  in  both 
systesM.  The  derivative  of  the  misali^nm.  t  sector  X  Is 


r  . 


(9) 
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The  vectors  Pm  and  p$  describe  the  angular  rates  of  master  and  clave 
platforms  with  respect  to  earth- fixed  coordinates  and  cm  and  es  are  the 
vectors  describing  respective  platform  angular  rater  due  to  gyro  drift. 

Equation  9,  resolved  in  the  master  coordinates,  can  be  expressed  as 


-(“s)£-(“m)m  +  c  *(“s), 

*(‘s  -  s)  +  (5s  -  «k)  +  «  •  (“s) 


(10) 


The  angular  rates  Pm  and  pg  are  determined  by  the  velocities  indicated 
in  the  inertia1,  systems,  i.e.. 
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Consequently, 


(13) 
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fepproriMtlni  as  an  angular  rate  vector  u  (roughly  equivalent  to 

fi  4-  P*  since  la  the  croaa  product  vita  the  small  vector  £  only  fire  - 
order  ter—  will  be  retained) ,  Eq .  10  and  3  3  can  be  coebined  to  give 

first— order  correct  expressions  for  the  components  of  £ 
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Equationr  7  and  14  can  be  written  in  vector-matrix  form  as 
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A  block  diagraa  Illustration  of  Eq.  15  Is  shown  in  Fig.  1. 
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The  measurements  provided  to  the  Kalman  Filter  are  the  velocity 
dlfferences  Avx,  Avy,  and  Avz.  The  measurement  matrix  for  the  state 
vector  in  Eq.  15  is 


(16) 


The  velocity  comparison  scheme  is  illustrated  in  Fig.  2. 


FIG.  2.  Velocity  Comparison  for  Kalman  Filter  Transfer  Alignment. 


Frequently  the  vertical  axis  of  near-earth  inertial  navigators  is 
instrumented  by  means  other  than  inertial  sensors.  This  is  a  consequence 
of  the  basic  instability  which  exists  in  computing  vertical  position  and 
velo'ity  from  accelerometer  measurements  (see  Appendix  C  of  Ref.  1). 

When  other  means  of  z-axis  instrumentation  are  provided,  Avz  becomes 
independent  of  the  other  quantities  in  the  state  vector  of  Eq.  15.  The 
third  row  in  the  matrix  contains  only  zeros.  The  third  column  in  the 
matrix  is  also  zero  and  Eq.  15  can  be  written  in  the  uncoupled  form 
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(18) 


Equation  18  is  of  no  value  in  determining  the  misalignment  angles  £x, 
iy,  and  E,z.  Equation  17  is  the  basis  for  the  linear  description  of 
error  (state)  dynamics  used  in  most  of  the  transfer  alignment  studies 
in  this  report.  The  corresponding  measurement  matrix  is  given  by 


H 


1  0  0  0  0 
0  10  0  0 


(19) 


The  gyro  and  accelerometer  errors  shown  in  Eq.  17  are  highly  cor¬ 
related  in  time.  Because  the  Kalman  Filter  formulation  presented  in 
Section  2  cf  Ref.  1  permits  only  uncorrelated  system  disturbances,  the 
state  vector  must  be  augmented  to  include  the  inertial  sensor  errors. 

It  is  assumed  that  the  time  to  perform  a  practical  transfer  alignment 
is  on  the  order  of  a  minute  or  two  and  the  gyro  and  accelerometer  errors 
will  remain  essentially  constant  over  that  period.  The  augmented  state 
vector  dynamics  are  described  by 
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where  I5  and  O5  are  the  5*5  identity  and  zero  matrices,  respectiveJ> , 
and 


V  -  V  -  7 
x  xs  "Si 


etc.  Notice  that  the  state  vector  differential  equation  (Eq.  20)  has  no 
forcing  terms.  The  corresponding  measurement  matrix  is 
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(21) 
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KALMAN  FILTER  ACCURACIES 

The  use  of  the  Kalman  Filter  to  estimate  the  state  variables  of  a 
linear  system  from  linear  measurements  permits  easy  calculation  of  the 
covariance  matrix  of  the  estimation  errors.  Because  system  velocity 
differences  are  provided  continuously,  the  matrix  Rlccati  differential 
equation  for  the  error  covariance  matrix  (Eq.  29  of  Ref.  1)  is  used 

P  -  FP  +  PFT  -  PHTR_1HP  +  GQ3T  (22) 


The  lack  of  random  forcing  terms  in  Eq.  20  permits  simplification  of 
Eq.  22 


P 


FP  +  PFT  -  PHTR-iHP 


(23) 


The  matrices  F  and  H  are  given  in  Eq  20  and  21.  The  matrix  R  and  initial 
conditions  on  the  error  covariance  remain  to  be  specified. 

Measurement  Errors 


The.  matrix  R  is  the  covariance  of  the  measurement  error  vector  v 
described  in  the  measurement  equation 


z  -  Hx  +  v 


(24) 


For  the  state  vector  and  measurements  described  under  Basic  Error  Dynamics 
(p.  2),  the  measurement  error  vector  v  is  composed  of  errors  in  indications 
of  the  differences  between  indicated  north  and  east  velocities  ox  the 
slave  and  master  inertial  platforms.  Velocity  differences  due  to  physical 
separation  of  tfie  two  platforms  are  remc.cd  by  compensating  the  measure¬ 
ment  under  the  assumption  that  the  connecting  structure  is  rigid.  Flexure 
of  this  structure  will  provide  uncompensated  random  velocity  differences, 
generally  of  a  high-frequency  nature  relative  to  the  system  response. 

These  errors,  in  addition  to  data  processing  and  signal  transmission 
noise  can  be  viewed  as  measurement  errors.  A  third  potential  velocity 
difference  error  arises  from  the  use  of  Doppler  velocity  measurements 
in  the  master  inertial  system.  All  three  effects  are  of  such  a  high- 
frequency  nature  relative  to  system  dynamics  that  they  will  be  considered 
uncorrelated.  Together  they  provide  the  measurement  error  v.  Because 
the  elements  of  v  are  north  and  east  velocity  difference  errors  while  the 
structural  vibration  and  Doppler  velocity  errors  are  usually  expressed  in 
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vehicle  axes ,  the  off-dtagonal  elements  of  R  are  generally  non-zero  and 
dependent  on  vehicle  orientation  relative  to  the  navigation  coordinates. 
However,  for  elmpllclty  the  two  elements  of  the  measurement  error  vector 
will  be  assumed  to  have  no  cross-correlation  in  the  study  which  follows. 
Under  the  assumptions  outlined  above,  the  combined  effects  of  vibration, 
data  noise,  end  Doppler  errors  will  be  described  by  a  single  quantity, 
the  diagonal  element  of  the  R  matrix,  r, 

fr  O'] 

R  -  (25) 

L° 

It  will  be  sufficient  for  the  purposes  of  this  study  to  vary  the  magnitude 
of  r  without  specifying  its  separate  contributing  sources.  (Indeed, 
measurement  error  sources  may  arise  in  addition  to  those  discussed  above.) 
For  example,  if  the  contributions  of  each  system  to  the  elements  of  v, 
vi,  and  V2  are  assumed  to  have  flat  spectral  density  with  bandwidth  of 
1,000  cps  and  RMS  magnitude  of  1.5  ft/sec,  Eq.  177  and  178  of  Ref.  1 
indicate  that 


r  a 


2(1.5  ft/sec)2  +  2(1.5  ft/sec)2 

(2 ir  x  io3  -M  (  2tt  x  103  — ) 
\  sec/  \  sec / 


m 


1.43  x  10 


(ft^/sec2)  sec 


Initialization 


Specification  of  the  initial  estimation  error  covariance  matrix  is 
simplified  by  the  assumption  that  there  is  no  cross-correlation  between 
estimation  errors  at  the  start  of  filtering.  This  is  a  reasonable  as¬ 
sumption  although  the  flexure  modes  of  the  carrying  vehicle  may  produce 
some  correlation  between  the  initial  values  of  misalignment  angles  £x, 

€y,  and  £z*  The  diagonal  elements  of  P(t0)  are  the  mean  square  errors 
in  the  initial  estimates  of  the  state  variables.  In  the  case  of  the 
velocity  differences  &vx  and  Avy,  it  is  assumed  that  the  master  velocities 
are  provided  for  initializing  the  slave  system  and  that  the  RMS  errors  in 
this  initialization  are  equal  to  the  RMS  errors  Introduced  in  the  velocity 
difference  between  the  two  systems.  For  example,  using  the  velocity 
indication  errors  described  above, 


^ll^o)  "  **22 (fco)  "  +  (1.5  ft/sec)^  •  4.50  ft2/sec2 
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The  mean  square  errors  In  the  initial  estimates  of  £x,  Kyt  and  £z 
are  taken  as  the  mean  square  misalignment  errors  between  the  two  systems. 
The  three  variances  are  chosen  to  be  equal  under  the  assumption  that  the 
factors  contributing  to  initial  misalignment  are  the  same  along  each  of 
the  navigation  system  axes.  For  example,  if  the  RMS  misalignment  errors 
were  two  degrees  about  each  axis. 


P33(£o)  ’  P«(£o)  ■  P55(£o)  ■  *  M'3 


The  diagonal  elements  of  P(t0)  corresponding  to  the  last  five  state 
variables  in  Gq.  20  can  be  easily  computed  from  the  mean  square  inertial 
sensor  errors.  Assuming  all  gyro  drift  rates  and  accelerometer  errors 
to  be  random  constants  with  zero  means  and  to  be  independent  from  one 
instrument  to  another. 


(26) 


where  V2g  and  are  the  mean  square  slave  x  accelerometer  and  master  x 
accelerometer  errors,  respectively.  The  remaining  diagonal  elements  are 
calculated  in  an  analogous  manner. 

Equation  26  illustrates  an  important  aspect  of  transfer  alignment. 

To  the  extent  that  Inertial  sensor  eriors  limit  transfer  alignment 
accuracy  or  practicality,  if  one  inertial  platform  contains  sensors  whose 
RMS  errors  are  much  larger  than  the  other  system's,  no  improvement  of  the 
more  accurate  set  of  inertial  sensors  will  improve  transfer  alignment 
significantly.  Accuracy  in  this  case  is  limited  by  the  quality  of  the 
poorest  set  of  sensors. 

As  an  example  of  the  calculation  of  the  lower  five  diagonal  elements 
in  the  initial  error  convariance  matrix,  consider  an  RMS  slave  accelerom¬ 
eter  bias  of  5  *  10"5  g  (1.61  *  1Q~3  ft/sec2)  and  an  RMS  bias  in  the 
corresponding  master  sensor  of  1C-'  j  (3.22  -  10-4  ft/sec?).  The  diagonal 
element  is  computed  according  to  Eq.  26 


-  ( 1.61  x  10"3)  +  1 3.  22  *  10"4)2  ss  ( 1. 61  x  l(f3)2 


_£  e  r  2 

-  2.59  x  10  — — T- 

4 

sec 

It  can  b®  seen  that  the  slave  sensor  error  dominates  these  terms. 
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The  eecCion  concludes  with  Che  observeclon  that  many  of  the  elements 
in  the  system  matrix,  F,  as  displayed  in  Eq.  20  may  be  time*- varying  and 
dependent  on  the  velocity  and  acceleration  of  the  vehicle  carrying  both 
Inertial  systems.  It  will  be  shown  later  that  these  elements  of  the  F 
matrix  play  an  important  role  in  determining  transfer  alignment  accuracy 
and  practicality. 
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Section  3.  IDENTIFICATION  OF  IMPORTANT  PARAMETERS 


A  large  number  of  parameters  exert  potential  Influence  on  the  Kalman 
Filter  accuracy  In  the  transfer  alignment  problem.  If  a  study  la  to 
proceed  efficiently  it  must  systematically  determine  which  parameters 
are  Important  over  the  range  of  values  considered  and  which  say  be  held 
constant  at  a  reasonable  value.  Once  the  important  parameters  have  been 
identified,  they  will  be  exercised  within  specified  limits  and  trade-off 
relations  will  be  determined. 


VELOCITY  AND  POSITION  MATCHING 

The  first  simplification  to  be  discussed  involves  transfer  alignment 
implementation  rather  than  inertial  system  parameters.  In  addition  to 
the  velocity-matching  scheme  discussed  in  Section  2,  it  is  also  possible 
to  compare  position  indications  from  the  two  inertial  systems  and  use  the 
differences  obtained  to  estimate  misalignment.  In  terms  of  the  matrices 
presented  in  Section  2,  this  scheme  adds  two  additional  rows  and  columns 
on  the  top  and  left  „,ide  of  F.  The  new  elements  are  all  zeros  except  for 
ones  in  the  third  column  of  the  first  new  row  and  the  fourth  column  of 
the  second  new  row.  The  H  matrix  also  contains  two  additional  rows  and 
columns.  They  are  added  on  the  right  and  lower  sides  and  consist  of  ones 
in  the  third  column  of  the  third  row  and  fourth  column  of  the  fourth  row, 
and  zeros  elsewhere.  The  R  matrix  is  expanded  to  dimension  4><4  and  the 
new  state  variables  Apx  and  Apy  are  placed  at  the  top  of  the  state  vector. 
Of  course,  the  error  covariance  matrix  also  contains  an  additional  two 
rows  and  two  columns. 

The  implementation  of  the  Kalman  Filter  in  the  case  where  position¬ 
matching  is  used  Involves  a  larger  F  matrix.  However,  the  significant 
Increase  in  complexity  results  from  doubling  the  number  of  Inputs  to  the 
filter  (measurements).  The  addition  of  two  new  state  variables  and  the 
new  measurements  raises  the  filter  gain  matrix  size  from  20  elements 
(2  *  10)  for  velocity-matching  to  48  elements  (4  *  12)  when  position- 
and  velocity-matching  are  used. 

To  evaluate  the  relative  accuracy  between  the  two  schemes,  the  error 
covariance  equation  (Eq.  23)  was  used  to  determine  RMS  alignment  accu¬ 
racies.  In  one  case  both  position  and  velocity  comparison  were  used, 
while  in  the  other  only  the  velocities  were  compared.  The  initial 
misalignments,  accelerometer  bias  errors,  and  velocity  measurement 
errors  used  are  those  presented  in  Section  2.^  The  RMS  slave  gyro  drift 


^Appendix  B  contains  a  table  of 
solution  of  Eq.  23  discussed  in  this 


Important  variables  for  each  computer 
document . 
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rat**  war*  taka*  to  be  0.05  deg/hr  and  the  aaiter  gyro  drift  rates  were 
assumed  much  nailer.  When  position-matching  was  used,  the  RhS  Initial 
error  in  position  difference  and  the  IMS  position  difference  measurement 
error  were  identical  to  the  values  for  the  corresponding  velocity  dif¬ 
ferences.  (Of  course  the  units  are  different.)  With  the  additional 
assumption  that  the  position  difference  errors  have  the  sare  frequency 
distribution  as  described  earlier  for  velocity  errors,  the  ^ laments  of 
&  and  P(0)  corresponding  to  position  differences  are  also  the  sane  as 
those  for  velocity  differences.  The  aircraft  carrying  the  two  systems 
was  assumed  to  be  in  level  flight  beading  due  east  with  a  velocity  of 
1,000  ft/sec.  (Aircraft  heading  has  no  Influence  on  alignment  accuracy. 

See  lef.  2.) 

Figures  3  and  4  show  the  SMS  estimation  error  histories  of  the  north 
level  and  azimuth  misalignment  angles,  respectively.  The  calculations 
represented  a  total  estimation  time  of  5  minutes  but  the  north  level 
error-  reached  its  lower  boundary  in  less  than  1  minute.2  E’st  level 
nlselignaent  behaved  in  the  same  manner.  The  combination  of  position- 
end  velocity-matching  provides  more  rapid  convergence  of  the  estimates 
on  the  true  valuee  but  does  not  result  in  greater  accuracy  than  velocity¬ 
matching  alone. 

While  the  accuracy  in  the  estimate  of  aslmuth  misalignment  Improves 
by  aura  than  an  order  of  magnitude,  the  SMS  value  of  £*  Is  still  much 
larger  than  those  for  the  two  level  misalignment  angles.  A  scheme  must 
be  sought  to  reduce  the  IMS  error  in  the  azimuth  misalignment  angle  to 
the  asms  magnitude  as  the  level  errors.  This  must  also  L*e  accomplished 
within  a  reasonable  time  period  if  1'  is  to  be  useful  frr  an  airborne 
missile.  In  the  next  section  it  will  be  demonstrated  thi*  greater 
accuracy  in  the  estimate  of  £z*  together  with  a  filter  time  on  the  order 
of  1  minute,  can  be  achieved  using  melocity-matching  alone.  Consequently, 
the  remainder  of  this  report  will  concentrate  on  this  much  simpler  scheme. 
The  combination  of  sufficient  accuracy,  rapid  convergence,  anc  relative 
slapliclty  make  the  velocity-matching  approach  considerably  nor  attractive 
than  the  comparison  of  both  position  and  velocity. 


BOklZOHTAL  MASEUVE2S 

If  Che  elrcreft  which  carries  the  two  Inertial  guidance  platforms 
is  permitted  to  perform  maneuvers  in  the  horizontal  plane  only,  the 
optimum  maneuver  is  a  high  acceleration  turn  (Kef.  2).  This  maneuver 
minimises  the  sum  of  the  mean  square  estimation  errors  of  the  three 
misalignment  angles.  The  major  benefit  of  the  turning  maneuver  Is  its 
ability  to  aid  the  Kalman  Filter  In  accurate  estimation  of  azimuth 


over  limit  of  the  RMS  north  level  estimation  error,  5  *  10~-> 
rad,  results  from  tbs  KSS  east  accelerometer  bias  error  of  5  *  10"^  g 
according  to  the  dlscussloo  In  Section  9  of  Bef.  1. 
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rIC.  3.  Coap*ri*on  cf  Position  and  Velocity  Matching  Vith 
Velocity  Matching  Kalman  Filter  Accuracies. 


3jH &**?;- 


me  TP  4653 


FIC.  4.  Comparison  of  Fosltiaa  aod  Velocity  Hatching 
ViCfc  Velocity  Matching  Kalun  Filter  Accuracies. 
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misalignment.  This  is  illustrated  in  Fig.  5.  As  the  horizontal  accel¬ 
eration  due  to  a  turning  maneuver  increases  from  0  to  3  gt  the  filter 
accuracy  increases  dramatically.  The  time  histories  shown  were  computed 
using  Eq.  23  with  the  initial  conditions  and  measurement  error  magnitude 
already  presented.  In  all  cases  the  turn  began  when  the  aircraft  was 
heading  due  south.  Reference  2  indicates  that  initial  heading  is  not 
significant.  Only  velocity  matching  was  considered.  Because  the  3 ~g 
horizontal  acceleration  (3 ,'l-g  total  acceleration)  is  considered  a 
reasonable  maximum  and  because  the  best  total  misalignment  accuracy  is 
obtained  with  that  maneuver,  all  results  reported  hereafter  will  be  for 
such  a  maneuver  unless  otherwise  stated.  In  addition,  all  parameters 
will  be  the  same  as  those  presented  above  unless  different  values  are 
given. 


GYRO  DRIFT  RATE 

The  effect  of  gyro  drift  rate  was  investigated  by  repeatfng  the 
calculations  for  the  y~g  turn  with  the  larger  RMS  slave  gyro  drift  rate 
of  0.5  deg/hr.  This  order-of-magnitude  larger  gyro  drift  rate  produced 
an  increase  in  the  sum  of  tne  RMS  misalignment  estimation  errors  at  the 
end  of  50  seconds  to  1.65  times  their  original  value.  Most  severely 
affected  was  the  RMS  azimuth  misalignment  error  which  rose  from 
1.04  *  10" 5  rad  to  3.12  *  10"^  rad.  However,  since  the  latter  val^c  is 
still  less  than  the  lower  limit  on  the  level  angle  estimation  errors, 
this  increase  is  less  significant  than  it  right  otherwise  seem-  The 
level  angle  estimation  errors  did  not  increase  as  dramatically  because 
they  had  essentially  reached  their  lower  limit  in  less  than  50  seconds 
when  the  original  gyro  drift  rate  ^alue  of  0.05  deg/hr  was  used.  Within 
the  range  described,  gy  ro  drift  rate  does  not  appear  to  exert  a  very 
significant  influence  on  the  Kalmar.  Filter  accuracy.  Consequently,  all 
subsequei r  results  will  be  obtained  assuming  the  lower  valu®  for  slave 
gyro  drift  at? . 


initial  misalignment  errors 

The  effect  of  the  initial  misalignment  errors  on  the  Kalman  Filter 
accuracy  and  filtering  time  was  studied  hy  increasing  the  RMS  initial 
misalignment  errors  (as  reflected  in  me  appropriate  diagonal  elements  of 
P(0))  from  2  to  4  deg/axis.  Again.  Eq.  23  vas  used  to  determine  tie" 
histories  of  error  statistics.  Nc  significant  Increase  was  found  in  the 
nif aligntaen'  estimation  errors  at  the  end  of  10  seronds.  It  is  concluded 
that  the  SMS  value  of  initial  misalignment  errors  has  no  effect  on  Kalman 
Til  ter  accuracy  in  the  parameter  range  of  interest.  The  RHS  magnit'ide  of 
2  deg'axis  will  be  used  in  all  subsequent  cal  col  at i ons . 
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FIG.  5.  Velocity  Matching:  Effect  of  Horizontal 
Acceleration  on  Kalman  Filter  Accuracy. 
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ACCELEROMETER  BIAS  ERROR 

As  already  discussed,  the  KMS  slave  accelerometer  bias  error  provides 
a  lower  limit  for  the  accuracy  of  Kalman  Filter  estimates  of  level  errors 
if  tha  aircraft  is  restricted  to  move  in  the  horizontal  plane.  Con¬ 
sequently,  it  ie  noc  surprising  that  raising  the  SMS  bias  value  by  an 
order  of  magnitude  increases  the  SMS  level  angle  estimation  errors  at 
the  end  of  50  seconds  by  a  factor  of  10.  Of  course,  the  simple  relation 
which  results  is  also  a  consequence  of  the  fact  that  in  both  cases  the 
errors  in  level  angle  estimates  have  reached  the  lower  limit  at  the  end 
of  50  seconds  of  filtering.  Figure  6  shows  the  north  level  error  time 
histories  for  the  two  bias  values. 3  No  significant  change  in  azimuth 
angle  estimation  accuracy  resulted  fiom  changing  the  level  accelerometer 
bias  errors.  Obviously,  at  least  for  the  case  where  the  aircraft  is 
restricted  to  horizontal  motion  during  alignment,  slave  accelerometer 
bias  errors  are  important. 


MEASUREMENT  noise 

It  is  of  particular  interest  to  observe  the  effect  of  variations  in 
measurement  noise  as  expressed  in  the  diagonal  elements  of  the  R  matrix. 
Not  only  is  the  measurement  noise  amplitude  to  be  encountered  in  a  prac¬ 
tical  application  not  well  known,  but  the  elements  of  R  represent  a 
combination  of  mean  square  measurement  errors  and  the  spectral  distribu¬ 
tion  of  these  errors.  Consequently,  variation  of  this  parameter  (diagonal 
element,  r)  is  necessary  to  provide  a  basis  for  predicting  filter  accuracy 
when  these  characteristics  are  subsenuenlly  determined.  Figures  7  and  8 
illustrate  the  effect  of  measurement  errors  on  the  ability  of  the  Kalman 
Filter  to  rapidly  provide  accurate  estimates  of  the  misalignment  angles. 
The  larger  value  of  the  R  matrix  diagonal  element,  r,  represents  measure¬ 
ment  errors  with  an  RMS  amplitude  8  times  that  represented  by  the  smaller 
value,  both  having  the  same  spectral  distribution.  The  larger  measurement 
error  causes  a  slower  convergence  of  the  Kalman  Filter  estimates  on  the 
true  values  of  the  misalignment  angles.  A  more  complete  presentation  of 
the  effects  of  measurement  errors  is  deferred  until  the  next  section. 


3 

Throughout  the  report  fjgures  illustrating  effects  on  level  angle 
estimation  accuracy  will  be  drawn  for  £x,  the  angle  about  the  north  axis. 
However,  all  conclusions  apply  also  to  £y  which  exhibits  almost  identical 
behavior. 
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VEHICLE  VELOCITY 

The  aircraft  forward  velocity  maintained  during  the  horizontal 
turning  maneuver  was  lowered  to  500  ft/aec — half  the  value  assumed  in  all 
other  computer  runs.  The  major  effect  of  changing  vehicle  velocity, 
assuming  the  same  horizontal  acceleration  is  maintained,  is  to  double 
the  frequency  of  the  sinusoidal  horizontal  acceleration  and  angular  rate 
components  in  the  F  matrix.  The  estimation  accuracies  computed  using 
Eq.  23  showed  only  a  slight  improvement  over  those  obtained  when  the 
velocity  was  1,000  ft/sec.  It  i.s  concluded  that,  over  this  range  of 
values,  the  forward  speed  of  the  aircraft  carrying  the  two  inertial 
platforms  have  little  effect  on  the  Kalman  Filter  alignment  accuracy. 


INSTRUMENT  ERROR  CALI BRAT ION 

When  the  Kalman  Filter  is  used  to  estimate  inertial  navigation  system 
errors,  the  possibility  exists  for  simultaneously  estimating  the  inertial 
sensor  errors — instrument  error  calibration.  Unfortunately,  during  all 
the  short-time  (less  than  »C  seconds)  alignment  runs  made  in  the  course 
of  this  investigation,  no  significant  improvement  in  the  accuracy  of 
sensor  error  estimates  was  obtained.  The  only  exception  during  the 
entire  study  occurred  when  position  and  velocity  differences  were  used 
for  a  period  of  5  minutes.  In  that  case  the  x-axis  and  y-axis  gyro  drift 
rates  were  reduced  by  a  considerable  factor.  The  possibility  of  Inertial 
sensor  calibration  during  Kalman  Filter  alignment  of  the  type  discussed 
appeals  remote.  Gyro  drift  rates  cause  errors  in  position  and  velocity 
to  grow  too  slowly  for  detection  in  a  reasonable  filtering  time  and  the 
initial  level  errors  are  so  large  relative  to  the  accelerometer  errors 
considered  that  no  improveo  estimate  can  be  made  when  a  horizontal 
maneuver  is  prescribed.  Section  5  discusses  the  possibility  of  using 
other  techniques  to  help  calibrate  sensor  errors. 

It  is  concluded  that  the  principal  parameters  which  determine  Kalman 
Filter  alignment  accuracy  when  the  aircraft  is  flying  in  a  horizontal  turn 
are  accelerometer  bias  errors  and  velocity  difference  measurement  errors. 
Consequently,  the  following  section  will  concentrate  on  establishing 
accuracy  trade-offs  involving  these  two  quantities. 
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Section  4.  TRADE-OFFS  INVOLVING  IMPORTANT  PARAMETERS 


ACCELEROMETER  BIAS  ERRORS 


The  accelerometer  bias  errors  were  varied  over  a  range  of  values 
between  5  *  10“ 5  and  5  *  10“*  g.  All  other  parameters  were  identical 
with  those  used  in  establishing  the  filter  accuracies  displayed  in  Fig.  5. 
The  aircraft  carrying  the  two  inertial  platforms  was  assumed  to  be  per¬ 
forming  a  horizontal  turning  maneuver  with  a  horizontal  acceleration  of 
3  g.  Figure  9  illustrates  the  RMS  north  level  estimation  error  for 
different  filtering  times  over  the  range  of  accelerometer  bias  errors. 

The  lower  limit  of  these  errors  as  a  function  of  bias  is  the  straight 
line  indicated  on  the  figure.  The  lines  representing  30,  40,  and 
50  seconds  of  filtering  reach  this  limit  as  the  bias  error  nears  the 
upper  bound  of  the  range  Investigated.  The  principal  effect  of  accel¬ 
erometer  bias  on  Kalman  Filter  accuracy  is,  of  course,  this  lower  bound 
on  the  level  estimation  accuracy  when  the  aircraft  is  restricted  to 
maneuver  in  the  horizontal  plane  (and  only  north  and  east  velocities  are 
compared).  (See  Section  5.)  The  fact  that  the  curves  representing  filter 
accuracy  for  different  filtering  times  slope  upward  indicates  that  larger 
accelerometer  bias  errors  increase  the  filtering  time  required  to  reach 
a  prescribed  accuracy.  As  mentioned  in  the  previous  section,  accelerom¬ 
eter  bias  errors  have  no  appreciable  effect  on  azimuth  angle  estimation 
accuracy . 


MEASUREMENT  ERRORS 

Measurement  error  amplitude  and  bandwidth,  as  reflected  in  the  R 
matrix  diagonal  elements,  affects  the  speed  with  which  the  Kalman  Filter 
estimates  converge  on  the  true  values  of  the  misalignment  angles.  A 
discussion  of  this  effect  is  found  on  pa*,’  20  of  Pef.  1.  Figures  10  and  11 
illustrate  the  accuracy  trade-offs  between  velocity  difference  measure¬ 
ments  end  the  Kalman  Filter  estimate  of  level  misalignment.  The  measure¬ 
ment  error  standard  deviation,  0V,  was  varied  between  1/2  and  8  times 
the  value  used  in  most  calculations  thus  far.  Again,  the  lower  limit 
imposed  on  estimation  accuracy  by  accelerometer  bias  errors  is  evident. 

For  e  bias  error  of  5  *  10~5  g  the  Kalman  Filter  accuracy  at  the  end 
of  10  seconds  is  linearly  related  to  the  RMS  measurement  noise,  ov. 

In  all  cases  the  differential  change  in  the  RMS  estimation  error  for  a 
given  filtering  time,  de,  is  related  to  a  differential  change  in  ov  by 


where  the  exponent  b  is  <1.  As  the  accelerometer  bias  error  or  the 
filtering  time  are  allowed  to  increase  within  the  range  investigated, 
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(ov  -  1.5  FT/SEC.  MEASUREMENT  NOISE  BAKDVIDTM  -  1,000  Hz) 


FIG.  9.  Velocity  Matching:  Effect  of  Accelerometer 
Bias  on  Kalman  Filter  Accuracy. 
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FTC.  10.  Velocity  Matching:  Effect  of  Measur  sent  Noise  on 
kalaan  Filter  Accuracy. 


HOC  TP  4653 


the  site  of  b  do'ruw*.  This  is  also  a  consequence  of  tbs  Halting 
offset  oa  Level  as t last lens  accuracy  provided  by  bias  errors:  If 
filtering  tine  and  accelerometer  error  are  such  that  the  limit  is  reached 
before  die  end  of  filtering,  a  considerable  increase  in  measurement 
errors  mast  take  place  h-’fore  the  level  errors  become  larger  than  the 
limit. 


Figure  12  illustrates  the  effect  of  varying  the  diagonal  elements 
of  the  1  matrix  on  azimuth  error  estimation  accuracy.  The  relations  are 
essentially  independent  of  accelerometer  bias.  The  curves  of  equal 
estimation  time  indicate  that  the  relation  between  azimuth  angle  estimation 
accuracy  aaa  IMS  measurement  errors  is  linear  for  fixed  filtering  duration 
and  measurement  error  bands ldth.  Since  the  parameter  varied  to  obtain 
these  relations  was  r ,  not  Oy,  the  accuracy  also  varies  inversely  as  the 
square  root  o-  beudwidth  when  ov  is  fixed.  Only  the  40-second  and 
50-seccr-4  curves  show  s  departure  from  the  first-order  relation.  They 
curve  at  the  low  values  of  indicating  a  somewhat  less  that  linear 
effect  In  that  region  of  parameters.  In  the  area  encompassing  small 
measurement  errors  and  long  filtering  times  the  accuracy  of  the  azimuth 
misalign— fit  estimate  aprroaches  a  non-zero  steady-stare  value  which  is 
not  dependent  on  measurement  noise  alone.  In  particular,  the  gyro  drift 
gates  are  also  important  in  establishing  the  steady-state  IMS  error  in 
it'  however,  over  almost  the  entire  range  of  av  and  filtering  time 
explored,  the  accuracy  of  the  azimuth  error  estimate  end  the  velocity 
difference  measurement  are  linearly  related  through  a  constant  determined 
by  the  measurement  error  bandwidth  and  the  duration  of  the  filtering 
process. 

Figures  9  through  12  provide  useful  trade-offs  for  preliminary 
determination  of  the  Kalman  Filter  alignment  accuracies  possible  with 

proposed  systems. 
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Section  5.  CONSIDERATION  OF  VERTICAL  MANEUVERS 


Depending  on  the  accuracy  desired,  the  Kalman  Filter  transfer 
alignment  of  a  slave  inertial  platform  with  parameters  within  the  ranges 
explored  in  this  report  may  not  be  satisfactory.  In  particular,  the 
lower  limit  or  level  error  estimation  accuracy  imposed  by  accelerometer 
bias  errors  may  be  restrictive.  Also,  depending  on  the  mission  time 
anticipated,  the  inability  to  calibrate  gyros  and  accelerometers  can  be 
prohibitive.  Techniques  exist  for  simultaneously  removing  most  of  these 
objections.  They  have  been  investigated  in  Ref.  2  and  found  to  be  very 
useful. 


VERTICAL  MANEUVERS 

The  use  of  a  vertical  maneuver  during  alignment  permits  the  separate 
observation  of  accelerometer  bias  and  level  error  in  the  horizontal 
velocity  differences.  Vertical  accelerations  generated  provide  a  time- 
varying  characteristic  for  the  (az  -  g)  element  in  the  system  matrix  of 
Eq.  17  where  formerly  only  the  constant  quantity,  g>  appeared.  The  level 
error  angles  then  generate  time-varying  contributions  to  the  measure¬ 
ments  Avx  and  Avy.  These  car*  be  distinguished  from  the  constant  con¬ 
tributions  of  accelerometer  bias  errors  without  increasing  the  size  of 
the  Kalman  Filter.  The  result,  as  described  in  Ref.  2,  is  the  removal 
of  the  lower  limit  on  level  error  estimation  accuracy  and  the  ability  to 
calibrate  the  accelerometers.  The  improvements  obtained  depend  to  a 
great  measure  on  the  relative  values  of  inertial  sensor  errors  and 
measurement  errors.  If  the  gyro  drift  rates  and  velocity  difference 
measurement  errors  permit  potential  filter  accuracies  much  better  than 
those  imposed  by  the  accelerometer  bias  during  3  horizontal  maneuver, 
considerable  improvement  results  from  providing  vertical  as  well  as 
horizontal  motion,  during  filtering. 

Figure  13  shows  the  reduction  in  the  s'.m  of  the  mean  square  errors  in 

A  A  A 

£x,  £y,  and  that  was  achieved  by  following  a  horizontal  turn  with  a 
sinusoidal  or  'porpoise"  maneuver  in  the  vertical  plane.  The  abscissa  is 
the  peak  vertical  acceleration  achieved,  armax.  At  the  start  of  the  ver¬ 
tical  man-'„  ’er  the  sum  cf  mean  square  errors  was  approximately  4  *  10~6 
rad^.  Most  of  the  reduction  was  achieved  in  the  estimation  errors  for 
level  angles  £x  and  £y.  The  addition  of  vertical  maneuvers  can  provide 
dramatic  ims  rovements  ir  alignment  accuracy. 
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z-AXIS  VELOCITY  COMPARISON 

If  vertical  accelerometers  are  available,  another  approach  exists 
for  removing  the  restriction  on  level  angle  estimation  accuracy.  Without 
employing  an  additional  vertical  maneuver,  the  z-axis  velocity  difference 
can  be  used  by  the  Kalman  Filter.  The  basic  formulation  is  based  on 
Eq.  15.  Time-varying  indications  of  level  misalignment  appear  separately 
in  the  vertical  velocity  difference.  Consequently,  they  can  be  separated 
from  the  contributions  of  accelerometer  bias  errors.  The  filter  gain 
matrix  of  the  complete  Kalman  Filter  contains  36  elements  instead  of  the 
20  needed  for  the  approach  used  in  Sections  3  and  4.  Increased  complexity 
is  traded  for  the  ability  to  achieve  more  accurate  misalignment  estimates 
during  a  horizontal  turn. 

Reference  2  indicates  that  the  results  of  z-axis  velocity  comparison 
are  essentially  the  same  as  for  the  vertical  maneuver — better  calibration 
of  accelerometers  and  removal  of  the  relation  between  bias  errors  and 
level  estimate  errors  discussed  in  Section  3. 
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Section  6.  CONCLUSIONS 


The  US'?  of  the  Kalman  Fi i  ter  to  permit  estimation  of  misalignment 
between,  two  inertial  guidance  platforms  has  been  explored.  Equations 
were  derived  describing  the  linear  dynamics  of  the  errors  involved  and 
all  the  necessary  matrix  a  were  developed. 

It  was  determined  that  horizontal  maneuvering  of  the  carrying  aircraft 
is  necessary  to  provide  accurate  estimates  of  azimuth  error  and  alignment 
times  on  the  order  of  a  minute.  The  most  important  parameters  within  the 
range  investigated  were  accelerometer  bias  errors  and  velocity  difference 
measurement  errors.  For  the  principal  scheme  explored,  the  accelerometer 
bias  imposes  a  lower  limit  on  the  filter's  ability  to  estimate  level 
misalignment  angles.  The  velocity  difference  errors  help  to  determine 
the  speed  with  which  the  filter  estimate  converges  on  the  true  values. 

In  particular,  a  linear  relation  was  found  to  exist  between  errors  in 
the  azimuth  angle  estimate  and  measurement  accuracies.  No  improvement 
in  inertial  sensor  calibration  resulted  when  the  basic  scheme  was  applied. 
Trade-off  curves  were  established  relating  the  misalignment  estimate 
accuracies  to  accelerometer  bias  errors,  measurement  errors,  and  filter 
times . 

Finally,  alternative  techniques  were  discussed  briefly  and  reference 
was  made  to  other  work  which  demonstrated  their  usefulness.  These 
approaches  involved  vertical  as  well  as  horizontal  aircraft  maneuvers 
and  the  comparison  of  z-axis  velocities  during  horizontal  turns.  It  is 
concluded  that  these  other  schemes  are  useful  but  somewhat  more  complex 
to  implement. 
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Appendix  A 

ALIGNMENT  ERRORS  IN  THE  MASTER  SYSTEM  AT  THE  TIME  OF 
TRANSFER  ALIGNMENT 

In  the  common  circumstance  where  the  slave  inertial  guidance  system 
is  of  poorer  quality  than  the  master  system,  all  the  errors  in  the  latter 
are  incorporated  into  the  former  at  the  conclusion  of  transfer  alignment. 
Unless  the  steering  commands  to  the  weapon  containing  the  slave  system 
are  always  provided  in  the  coordinate  frame  defined  by  the  master  sy»tem 
(and  therefore  Include  errors  prcvided  to  the  slave  during  alignment) , 
the  nature  of  master  system  errors  at  the  time  of  transfer  alignment  is 
Important.  For  this  reason  one  set  of  such  errors  was  computed  for  the 
situation  in  which  the  aircraft  had  been  cruising  for  1  hour  prior  to 
initiation  of  the  alignment  procedure.  It  is  emphasized  that  this  set 
of  master  system  errors  is  highly  dependent  on  the  inertial  sensor 
errors  postulated,  and  especially  on  the  kinds  and  method  of  combination 
of  external  navigat  on  information  used  to  improve  the  master  Inertial 
system  accuracy  during  flight.  The  only  aid  assumed  in  this  example  is 
Doppler  velocity  information.  This  external  measurement  is  eimloyed  in 
an  optimal  manner  through  the  use  of  a  Kalman  Filter.  Though  in  actual 
practice  the  Doppler  information  may  be  time-averaged  before  use  (Ref.  3) , 
for  the  calculations  presented  here  it  is  considered  as  a  continuous 
measurement  with  high-frequency  noise. 

The  carrying  aircraft  is  assumed  to  have  cruised  due  east  in  level 
flight  with  c  velocity  of  1,000  ft/sec  for  1  hour  prior  to  alignment. 

The  following  parameter  descriptions  were  used: 

Latitude,  deg  .  45 

Doppler  measurement  noise  (unbiased),  ft/sec  RMS  .  3 

(across-  and  along-track  with  a  bandwidth  of  1  Hz) 

Initial  master  misalignment  angles,  mrad  RMS  . .  2 

Initial  master  position  error,  ft  RMS  .  10 

Initial  master  velocity  error,  ft/sec  RMS  .  2 

Master  accelerometer  error,  g  RMS  .  10-^ 

(exponentinlly-correlated  with  a  10-hr 
correlation  time) 

Master  gyro  drift  rate,  d  ^/br  RMS  . .  10~^ 

(exponentially-correlated  with  a  10-hr 
correlation  time) 
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The  error  covariance  matrix  differential  equation  was  solved  by 
numerical  Integration  and  at  the  end  of  1  hour  of  flight  the  RMS  errors 
in  the  master  inertial  navigator  were  calculated  to  be: 

Position,  ft 

North  . 

East  . 

Velocity,  ft/sec 

North  . 

East  . 

Level,  rad 

North  . .  5.0  x  10”^ 

East  . . .  1.1  x  10“  5 

Azimuth,  rad . . .  1.1  x  10-^ 

The  large  position  errors  occur  because  no  external  position  fixes  were 
made  during  the  flight.  If  initialization  of  the  slave  system  with 
respect  to  position  in  earth-fixed  coordinates  is  important,  some  con¬ 
sideration  must  be  made  for  providing  external  position  measurements  to 
the  master  navigator  prior  to  transfer  alignment.  Of  course,  addition 
of  position  measurements  would  only  serve  to  further  reduce  the  other 
master  system  errors  as  well.  No  significant  inertial  sensor  error 
calibration  occurred  during  this  ^.ruise  flight. 

The  fact  that  RMS  errors  in  all  other  quantities  decreased  during 
the  course  of  the  1-hour  flight  illustrates  the  self-alignment  capability 
of  the  Doppler  velocity  aided  master  inertial  navigator.  It  can  be  seen 
that,  for  the  range  of  values  investigated  in  Section  4,  the  transfer 
alignment  attitude  errors  are  of  the  same  order  of  magnitude  as  those 
already  present  in  the  master  inertial  system  at  the  start  of  alignment. 
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Appendix  B 
TABLE  OF  PARAMETERS 


Table  1  shews  the  values  of  all  important  parameters  specified  in 
the  computer  solutions  of  Eq.  23  which  are  discussed  in  this  report. 

The  section,  and  in  most  cases  the  figure,  in  which  the  results  for  each 
parameter  group  are  summarized  are  listed  in  the  first  column. 


TABLE  1.  Parameter  Groups 
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